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We present the first ab initio coupled-cluster calculations of medium-mass nuclei with explicit chiral three- 
nucleon (3N) interactions. Using a spherical formulation of coupled cluster with singles and doubles excitations 
including explicit 3N contributions, we study ground states of 16 ' 24 0, 40 - 48 Ca and 56 Ni. We employ chiral NN 
plus 3N interactions softened through a similarity renormalization group (SRG) transformation at the three-body 
level. We investigate the impact of all truncations and quantify the resulting uncertainties — this includes the 
contributions from triples excitations, the truncation of the set of three-body matrix elements, and the omission 
of SRG-induced four-body interactions. Furthermore, we assess the quality of a normal-ordering approximation 
of the 3N interaction beyond light nuclei. Our study points towards the predictive power of chiral Hamiltonians 
in the medium-mass regime. 
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The ab initio description of medium-mass nuclei is one of 
the most dynamic frontiers in nuclear structure theory today — 
bridging the gap between accurate ab initio calculations for 
light nuclei and the realm of approximate or phenomenologi- 
cal approaches for heavy nuclei and nuclear matter. A number 
of many-body methods are being developed and extended to- 
wards the medium-mass regime. Coupled-cluster (CC) theory 
has a pioneering role in this domain and other meth- 

ods, like self-consistent Green's function approaches J3l or 
the in-medium similarity renormalization group JHBl are fol- 
lowing. Extensions of the no-core shell model (NCSM) J3], 
like the importance-truncated NCSM HI @], are connecting 
the domains of light and medium-mass nuclei. 

A critical ingredient for all ab initio many -body approaches 
is the Hamiltonian. At present, chiral effective field the- 
ory (EFT) provides the most systematic approach to QCD- 
based Hamiltonians for accurate nuclear structure calculations 
[Till . Already the present generation of chiral Hamil- 
tonians, consisting of nucleon-nucleon (NN) interactions at 
next-to-next-to-next-to-leading order (N 3 LO) II12[ 1 1 3fl and 
three-nucleon (3N) interactions at next-to-next-to-leading or- 
der (N 2 LO) |l4|l give a very good description of p-shell nuclei 
as demonstrated in ab initio NCSM calculations 

EH. On- 
going developments in this sector, e.g., regarding chiral 3N 
interactions at N 3 LO IU9U21I1 or a A-full formulation of chiral 
EFT 112211 . will soon provide next-generation chiral Hamilto- 
nians with consistent NN and 3N interactions. 

The inclusion of the 3N interaction is vital to realize the pre- 
dictive potential of chiral Hamiltonians, but poses a number of 
computational challenges. So far, there are no calculations for 
medium-mass nuclei that include explicit chiral 3N interac- 
tions, without resorting to approximate or even schematic re- 
ductions to effective two-body interactions J3, 23]. A system- 
atic though approximate inclusion of chiral 3N interactions for 
medium-mass nuclei is the normal-ordering scheme applied in 
Ref. J24|. 

In this Communication we present, for the first time, ab 
initio coupled-cluster calculations for medium-mass nuclei 



including explicit 3N interactions. We have developed a 
spherical implementation of coupled cluster with singles and 
doubles excitations for three-body Hamiltonians (CCSD3B) 
which enables us to perform converged ground-state calcu- 
lations for closed-shell nuclei with full three-body interac- 
tions. In this framework we study the ground-state ener- 
gies of 16 ' 24 0, 40 - 48 Ca and 56 Ni using chiral NN+3N Hamil- 
tonians softened throu gh similarity renormalization group 
(SRG) transformations 1 1 8l 12511 . We systematically address 
all truncations introduced in the many-body framework and 
the Hamiltonian and quantify the resulting uncertainties in 
the ground-state energies. We demonstrate that the overall 
uncertainty for the prediction of ground-state energies in the 
medium-mass regime is of the order of a few percent. Within 
these uncertainties the chiral NN+3N Hamiltonians used in 
this work predict ground-state energies that are in agreement 
with experiment. This is a remarkable result, since no infor- 
mation beyond the few-body domain (A < 4) was used to fix 
the parameters of the chiral interactions. 

Coupled-Cluster Method. In single-reference CCSD the 
ground state |*P) of a many-body Hamiltonian is parametrized 
by the exponential ansatz 
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with T = T\ + T2, where T„ are excitation operators of the 
form 
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acting on a single Slater-determinant reference state |<I>). 

Coupled-cluster theory is conveniently formulated in terms 
of Hamiltonians that are normal ordered with respect to |<J>). 
If the original Hamiltonian 



H - hi + h 2 + hi 
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consisting of one-, two- and three-body contributions is cast 
into normal-ordered form, 



H = <<D|//|4>> + h™° + t%° + h^° 
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contributions from the n-body part h n of the original Hamil- 
tonian enter the matrix elements of the normal-ordered opera- 
tors of particle ranks k <n. 

Neglecting the three-body part of the normal-ordered 
Hamiltonian (O and thus only including contributions of the 
three-body interaction that have been demoted to lower parti- 
cle ranks through normal-ordering leads to the N02B approx- 



spherical scheme is given by 



imation discussed in 1241 12611 . In this Communication we in- 



clude the complete 3N interaction, i.e., we include the normal- 
ordered three-body part beyond the N02B approximation 
explicitly. The CCSD energy and amplitude equations includ- 
ing the full 3N interaction can be written as 
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where A£no2b and T„ ,no2b denote the standard CCSD equa- 
tions for two-body Hamiltonians i27ll corresponding to the 
N02B approximation. Thus, the inclusion of the residual 
three-body operator /z? generates a set of additional terms 
in the CCSD equations. 

Expressions for (E])-© in ;7?-scheme are presented in a 
factorized form in 12611 . We rederived these equations in a 
straightforward unfactorized way, resulting in more but sim- 
pler terms. Already for two-body Hamiltonians the basis 
sizes and particle numbers for m-scheme CC calculations are 
severely limited by the number of amplitudes and matrix el- 
ements that need to be handled. It is well known that the 
range of the CC method can be greatly extended by exploit- 
ing spherical symmetry and using an angular-momentum cou- 
pled formulation 112 811 . We have developed such an efficient 
spherical implementation of CCSD3B. Proceeding along the 
lines of Refs. 12811 . we couple the external lines of the dia- 
grams, cut open internal lines and perform angular momen- 
tum recouplings in order to express diagrams in terms of 
angular-momentum coupled matrix elements of the operators 
involved. For example, one of the computationally more in- 
volved contributions to the T2 amplitude equations reads in 
m-scheme 
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where we use the standard diagrammatic representation 112711 
and the permutation operator P pq — 1 - T pq with T pq denoting 
index transpositions. The corresponding contribution in the 
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T pq and x 



— ^l2x + 1. Here, indices 
a,b, . . . represent («Zy')-orbitals rather than individual single- 
particle states. The coupling lines indicate standard angular- 
momentum coupling and the tilde denotes time-reversed or- 
bitals. 

In our spherical scheme, the three-body matrix elements of 
enter in a specific reduced and coupled form, given by 

3 ah J(ie 

(S>c\\hT\&f) = X' ' ' J " j: { 'i 1i j j } 

I i J 

x ([(ab)J ab c]JM\h™°\[(de)J de f]JM) , (10) 

where the 3N matrix element on the r.h.s is M-independent. 
We use this specific coupling to store the three-body matrix 
elements in fast memory, making use of a variant of the ef- 
ficient storage and retrieval schemes we have developed for 
,/T-coupled three-body matrix elements 1 1 8l 12911 . This is crit- 
ical for the over-all performance of the CCSD3B calculations. 
In order to accelerate the convergence of the iterative solution 
of the CCSD3B amplitude equations, we initialize the ampli- 
tudes with the solution of the corresponding CCSD calcula- 
tion in N02B approximation. In addition to the CCSD and 
CCSD3B calculations we will perform A-CCSD(T) calcula- 
tions lEil [30I l3lll using the N02B Hamiltonian to assess the 
influence of triples excitations. 

Hamiltonian and Basis. The starting point for our in- 
vestigation of medium-mass nuclei are SRG-evolved chiral 
NN+3N Hamiltonians. We use the chiral NN interaction at 
N 3 LO |l2| and a local form of the chiral 3N interaction at 
N 2 LO UM. Instead of a momentum cutoff of 500 MeV used, 
e.g., in Ref. 118I1 . we reduce the cutoff of the initial 3N interac- 
tion to 400 MeV and choose eg = 0.098 to reproduce the 4 He 
ground-state energy, keeping cd = -0.2 Il24ll . This cutoff re- 
duction is motivated by the observation that SRG-induced 4N 
interactions have a sizable impact on ground-state energies of 
medium-mass nuclei, which can be reduced efficiently by re- 
ducing the cutoff of the initial 3N interaction jl8Ll24l29ll . We 
emphasize that the following results for medium-mass nuclei 
from ls O to 56 Ni are pure predictions. 

We will employ two types of SRG-evolved Hamiltonians: 
The NN+3N-full Hamiltonian starts with the initial chiral 
NN+3N Hamiltonian and retains all terms up to the 3N level 
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in the SRG evolution, the NN+3N-induced Hamiltonian omits 
the chiral 3N interaction from the initial Hamiltonian, but 
keeps all induced 3N terms throughout the evolution. In ad- 
dition, we use a range of flow parameters a in order to assess 
the role of SRG-induced contributions beyond the three-body 
level d]. 

The underlying single-particle basis is a harmonic- 
oscillator basis truncated in the principal oscillator quantum 
number 2« + / = e < e max and we go up to e max = 12. We 
perform Hartree-Fock (HF) calculations including the 3N in- 
teraction for each set of basis parameters to obtain an opti- 
mized single-particle basis and stabilize the convergence of 
the CC iterations. The normal-ordering is done consistently, 
i.e., with respect to the HF reference state. At the moment 
it is not possible to include all three-body matrix elements 
that would appear in the larger bases, we are limited to three- 
body matrix elements with e\ + e2 + < Ej, m - dx = 14. We 
will discuss the impact of this additional cut in detail later on. 
We stress that the exact treatment of the isospin dependence 
of the two- and three -body matrix elements in the HF basis 
is crucial. It is generated by isospin-dependence of the HF 
single-particle wavefunctions, although the 3N operator used 
here is isospin symmetric. In order to avoid a drastic increase 
of storage needed for the HF three-body matrix elements, we 
perform the transformation to the HF basis on the fly. 

Results. We present first results of CC calculations with 
explicit 3N interactions for the ground-state energies of 16,24 0, 
40 ' 48 Ca, and 56 Ni. Step by step, we will quantify the uncertain- 
ties resulting from various truncations in the many-body treat- 
ment. First we contrast CCSD3B calculations with CCSD us- 
ing the N02B approximation for the Hamiltonian. For both 
we discuss basis-space convergence in terms of e max and os- 
cillator frequency tiQ.. We then study the influence of the 
three-body truncation E^ m!ix . Finally, we include non-iterative 
triples corrections at the level of A-CCSD(T) with the N02B 
Hamiltonian. In all cases we vary the flow parameter a over a 
wide range to study the impact of induced many-body terms. 

Figure Q] shows a comparison of CCSD3B calculations us- 
ing the complete 3N interaction with CCSD using the N02B 
approximation for the ground-state energies of 16,24 0, 40 ' 48 Ca, 
and 56 Ni as function of the basis size e max for, both, the 
NN+3N-induced and the NN+3N-full Hamiltonian with three 
different values of the SRG flow parameter. The oscillator fre- 
quencies correspond to the energy minima in the largest basis 
spaces (cf. Fig. 0). The numerical values of the ground-state 
energies for the largest basis sets are also summarized in Tabs. 
H] and |n] The first observation is that we are able to converge 
or come sufficiently close to convergence with respect to the 
basis size e max in practically all cases. The second observation 
is that the N02B works extremely well for all cases: For 16 
the largest deviation from the full CCSD3B results is 0.9 Me V 
or 0.8%, for 56 Ni it is 4 MeV or 0.8% across all Hamiltonians 
considered here. Given that the computational cost for the 
CCSD3B calculations is two orders of magnitude higher than 
for CCSD and that the accuracy we target for many-body cal- 
culations in this mass range is not better than 1%, the N02B 



approximation constitutes a very efficient tool. 

The quality of the N02B approximation is confirmed in 
Fig. |2] where it is compared to CCSD3B using the NN+3N- 
full Hamiltonian and the largest basis set as function of the os- 
cillator frequency. The accuracy of the N02B approximation 
is largely independent of TiQ.. Note that the effect of the resid- 
ual 3N interaction beyond the N02B approximation is always 
repulsive, i.e., of the same sign as the complete 3N contribu- 
tion composed of induced and evolved initial 3N terms 12411 . 

The fact that the SRG evolution in the 3N sector is per- 
formed in a finite model space of harmonic oscillator Jacobi 
states |18, 2^] leads to additional uncertainties at low frequen- 
cies fiQ.. By varying the size of SRG model space 12911 we 
estimate the uncertainties at the optimal frequency to be much 
smaller than 1 % for all nuclei except 56 Ni, where they reach 
the 1 % level. For smaller frequencies, however, the truncation 
of the SRG model space leads to more significant effects — the 
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FIG. 1: (color online) Ground-state energies for 1624 0, 40,48 Ca and 
56 Ni as function of e raax for the two types of Hamiltonians (see col- 
umn headings) using CCSD3B (solid lines) and the N02B approxi- 
mation (dashed lines) for a range of flow parameters: a = 0.02 fm 4 
(•), 0.04 fm 4 (♦), 0.08 fm 4 (a). 
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TABLE I: Summary of CCSD and A-CCSD(T) ground-state energies 
in MeV for the NN+3N-induced Hamiltonian for a subset of a values 
computed at optimum oscillator frequency tiCl = 20 MeV for 16 ' 24 0, 
tiQ. = 24 MeV for 40 Ca, and hQ, = 28 MeV for 48 Ca, 56 Ni. The 
number in parenthesis denotes the i?3 max cut in the 3N Hamiltonian. 
The last column gives e max -extrapolated values, the other columns 
are for e max = 12. 



NN+3N 




CCSD 


CCSD 


A-CCSD(T) 


induced 


o-[fm 4 ] 


3B 


N02B 


N02B 


N02B 


N02B™ 






(12) 


(12) 


(14) 


(14) 


(14) 


16q 


0.02 


-121.6 


-122.3 


-121.7 


-126.1 


-126.5 




0.04 


-121.1 


-121.9 


-121.6 


-124.4 


-124.4 




0.08 


-119.9 


-120.8 


-120.8 


-122.4 


-122.4 


24 Q 


0.02 


-156.1 


-157.0 


-155.8 


-162.7 


-163.7 




0.04 


-155.4 


-156.6 


-155.9 


-160.2 


-160.4 




0.08 


-153.3 


-154.5 


-154.4 


-156.9 


-157.1 


40 Ca 


0.02 


-362.1 


-363.9 


-360.5 


-374.3 


-375.5 




0.04 


-357.7 


-359.9 


-358.1 


-366.3 


-366.6 




0.08 


-348.8 


-350.9 


-350.8 


-355.3 


-355.5 


48 Ca 


0.02 


-428.1 


-430.5 


-425.3 


-442.3 


-443.9 




0.04 


-419.6 


-422.3 


-420.0 


-429.6 


-430.9 




0.08 


-403.4 


-405.7 


-406.8 


-411.9 


-413.3 


56 Ni 


0.02 


-473.2 


-475.8 


-464.3 


-487.5 


-489.8 




0.04 


-462.6 


-465.6 


-458.0 


-472.4 


-474.2 




0.08 


-439.3 


-441.8 


-439.3 


-448.4 


-450.4 


TABLE II: Same as in Tab.Ufor the NN+3N-full Hamiltonian. 


NN+3N 




CCSD 


CCSD 


A-CCSD(T) 


full 


a[fm 4 ] 


3B 


N02B 


N02B 


N02B 


N02B 00 






(12) 


(12) 


(14) 


(14) 


(14) 


16q 


0.02 


-124.9 


-125.4 


-124.5 


-129.9 


-130.4 




0.04 


-126.8 


-127.6 


-127.1 


-130.8 


-130.8 




0.08 


-128.2 


-129.1 


-129.0 


-131.2 


-131.2 


24 Q 


0.02 


-161.4 


-162.0 


-160.4 


-168.3 


-169.4 




0.04 


-165.3 


-166.4 


-165.5 


-170.7 


-170.9 




0.08 


-167.6 


-168.8 


-168.6 


-171.6 


-171.7 


40 Ca 


0.02 


-354.4 


-356.0 


-352.1 


-370.6 


-371.7 




0.04 


-363.8 


-366.2 


-364.3 


-376.4 


-376.7 




0.08 


-369.0 


-371.5 


-371.3 


-378.2 


-378.4 


48 Ca 


0.02 


-429.4 


-431.6 


-426.7 


-449.5 


-450.9 




0.04 


-441.2 


-444.3 


-441.9 


-456.0 


-457.0 




0.08 


-445.3 


-448.2 


-448.3 


-453.5 


-456.7 


56 Ni 


0.02 


-497.3 


-499.9 


-490.9 


-521.7 


-523.4 




0.04 


-513.9 


-517.9 


-511.4 


-530.9 


-531.8 




0.08 


-517.0 


-520.7 


-517.9 


-528.4 


-529.2 



increase of the ground-state energies of Ca and beyond at 
the lowest frequencies shown in Fig. |2]is partly due to this. 

Next we address the £3 max cut used in the 3N matrix ele- 
ments for technical reasons. In Tabs. |T]and|ll]the CCSD results 
using the N02B approximation with £3„ lax = 12 and 14 are 
compared. We find that the influence of this cut grows with 
increasing particle number and decreasing flow parameter a. 
For the softest interaction with a = 0.08 fm 4 the cut is com- 
pletely irrelevant up to 40 Ca, only for 56 Ni we observe a 0.5% 
change in the ground-state energy. For less evolved interac- 
tions the effect increases and reaches about 1 % up to 48 Ca and 



about 2% for Ni. For the description of still heavier nuclei 
or the use of bare 3N interactions one will have to improve on 
this truncation in order to reach accurate results. 

As the final model-space related truncation, we discuss the 
truncation of the excitation operators in the coupled cluster 
ansatz (Q]). In Fig. [3] we compare the results of CCSD and 
A-CCSD(T) calculations for all nuclei and Hamiltonians in 
N02B approximation as function of the basis truncation e max 
for £3 m ax = 14. The results for e max = 12 are also summa- 
rized in Tabs. |l]and|ll] Again we observe a systematic depen- 
dence on the flow parameter a. For the softest interactions 
with a = 0.08 fm 4 the inclusion of triples excitations lowers 
the ground-state energies by 1.5 to 2% for all nuclei and both 
Hamiltonians. For a = 0.02 fm 4 the difference increases to 
about 4 to 6%. If we conservatively consider the triples con- 
tribution as a measure for the inherent uncertainty due to the 
truncation of the cluster operator, then this is the largest un- 
certainty so far. 

Finally, we quantify the uncertainty due to the omis- 
sion of the SRG-induced four- and more-nucleon interactions 
through the a-dependence of the CCSD and A-CCSD(T) re- 
sults shown in Fig. [3] and Tabs. |T] and [II] First of all, we 
note that missing many-body terms of the Hamiltonian are 
of opposite sign but of the same order of magnitude for the 
NN+3N-induced and the NN+3N-full Hamiltonian. For the 
NN+3N-induced Hamiltonian the ground-state energies for 
a = 0.02 fm 4 (harder interaction) are systematically lower 
than for a = 0.08 fm 4 (softer interaction). Furthermore, the 
energy spread over this range of flow parameters is smaller for 
CCSD and larger for A-CCSD(T). For 40 Ca, e.g., the spread 
amounts to 3% of the ground-state energy in CCSD and 5% 
in A-CCSD(T). The pattern is reversed for the NN+3N-full 
Hamiltonian. The ground-state energies for or = 0.02 fm 4 are 
systematically above the energies for a = 0.08 fm 4 and the en- 
ergy spread is reduced by including the triples correction. For 
40 Ca the relative energy spread is 5% for the CCSD and 2% for 




20 24 28 32 36 20 24 28 32 36 
fiil [MeV] m [MeV] 



FIG. 2: (color online) Ground-state energies for 24 , 40,48 Ca and 
56 Ni as function of TiQ. at e max = 12 for the NN+3N-full Hamiltonian 
using CCSD3B (solid lines) and the N02B approximation (dashed 
lines) for a range of flow parameters: a = 0.02 fm 4 (•), 0.04 fm 4 ( 
♦), 0.08 fm 4 (A). 
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A-CCSD(T). Note that these a-dependencies are distorted by 
the influence of the Z^max-truncation. With increasing £3 max 
the A-CCSD(T) ground-state energies will move up for harder 
interactions and, thus, the apparent ff-dependence will be re- 
duced for NN+3N-induced and increased for NN+3N-full. 

If we compare the ground-state energies throughout the set 
of nuclei discussed here with experiment (cf. Fig. [3]), keeping 
in mind the uncertainties we discussed above, then the agree- 
ment is remarkable. We stress that no information beyond 
A = 4 was used to constrain the Hamiltonian, so obtaining 
the correct binding systematics for medium-mass nuclei is far 
from trivial. Though the impact of the initial chiral 3N interac- 
tion is moderate on the scales shown in Fig. [3j it is important 
to obtain the correction binding-energy systematics along iso- 
topic chains. In contrast, the effect of the SRG-induced 3N 
interactions is huge and their inclusion is mandatory — for 
56 Ni and a — 0.08 fm 4 a calculation with only SRG-evolved 

NN+3N induced NN+3N full 




4 6 8 10 12 4 6 8 10 12 



^max ^max 

FIG. 3: (color online) A-CCSD(T) (solid lines) and CCSD (dashed 
lines) ground-state energies for 16 - 24 0, 4048 Ca and 56 Ni as function of 
e m ax for the two types of Hamiltonians (see column headings) using 
the N02B approximation for flow parameters: a = 0.02 fm 4 (•), 
0.04 fm 4 (♦), 0.08 fm 4 (a). 



NN interactions yields an unphysical ground-state energy of 
about -950 MeV. 

Conclusions. We have presented the first ab initio calcu- 
lations for nuclei in the medium-mass regime with inclusion 
of explicit 3N interactions. On this baseline we quantify the 
effects of all truncations in the many-body approach, i.e., the 
cluster-rank of the CC ansatz, the single-particle truncation 
e m ax, the truncation of the 3N matrix elements £3 max, and the 
optional omission of residual normal-ordered 3N terms. For 
all truncations we clearly benefit from the prediagonalization 
of the Hamiltonian through the SRG evolution — for Hamil- 
tonians with a = 0.08 fm 4 the uncertainties due to the £3 max 
truncation and the N02B approximation are at the level of 
1% or below and the effect of the triples correction is only 
2%. In particular the impact of the £3 ma x truncation and the 
triples corrections increase rapidly when going to harder inter- 
actions. The omission of SRG-induced four- and many-body 
terms, i.e. the truncation of the particle rank of the Hamil- 
tonian, introduces an uncertainty at the level of a few % and 
thus limits the over-all accuracy of the approach. A reduction 
or elimination of this uncertainty will be a prime goal of future 
studies. Finally, our results point towards the predictive power 
of chiral NN+3N Hamiltonians for medium-mass nuclei. 
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